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We discuss a new mechanism which could cause a rotation of polarization of electromagnetic 
waves due to magnetic fields on cosmological scales. The effect is due to the geometrical phase of 
Pancharatnam and Berry, and causes a corkscrew twisting of the plane of polarization. The new 
effect represents an additional tool that allows possible intergalactic and cosmological magnetic fields 
to be studied using radio propagation. 

PACS numbers: 98.80.Es, 41.20.Jb 

Many studies have been done to find a cosmological component to magnetic fields. A common method looks for sys- 
tematic dependence of the Faraday rotation measures (RM) with increasing redshift z. For a radio wave of wavelength 
A and polarization angle Faraday rotation [jD predicts i9 = RMA 2 + %■> where RM = [e 3 /(167r 3 TO e £oc 3 )] J z B z n e dz 
for propagation through a magnetic field with component B z along the propagation path; n e is the electron number 
density of the plasma, and \ is t ne angle of emission. 

Observation of the Faraday effect measures integrated magnetic plasma parameters along the line of sight. Although 
extremely stringent bounds have occasionally been stated Q, the methodology has its limitations, inasmuch as reversal 
of the magnetic field will reverse the rotation measure. This is not always discussed. When possible magnetic field 
reversals are taken into account, the observed rotation measures are not capable of bounding local magnitudes of the 
fields: magnetic fields larger than some bounds could well exist if the fields go through some reversals of sign along 
the way. 

Rotation of the plane of polarization is usually attributed to birefringence. Here, we present a new physical 
phenomenon, which we call the "corkscrew effect," which could also create a rotation of polarization. The underlying 
observation is that any phase shift of the positive relative to the negative helicity amplitudes of the wave will rotate a 
linear polarization. The corkscrew effect occurs for propagation with angular frequency uj in an electromagnetic wave 
equation of the form 

[tfS^-kikj-^eijiuj^Ej =0; e^uj) = egM + (!) 

where e s (cu) is real and symmetric, and e A (oj) is non-zero, real and antisymmetric; the full dielectric constant is e(uj) 
Hermitian (no absorption) || . We are interested in circumstances where the propagating transverse modes are very 
nearly degenerate (free space Sij = Sij in our units), with small perturbations lifting the degeneracy. The system then 
becomes a two-state evolution problem, analogous to two-state mixing problems encountered in quantum mechanics, 
with e(w) with e representing parameters that can vary slowly along the path of propagation. Where slow parameter 
variations are an issue, naive inferences about wave propagation can be wrong. For example, if k\(z) and £2(2) are 
the wave numbers from diagonalizing the dispersion relation for fixed parameters and fixed to, one does generally not 
find the phase of a wave from the expected eikonal factors J i{k\(z), k2{z)} dz. This is due to the geometric phase 
phenomenon from parallel transport in the adiabatic limit, as described for visible-light optics by Pancharatnam and 
Berry Q. We focus here on this purely geometrical effect, noting for later study that strong dynamical effects can 
also occur - such as seen in resonant magnetic moment transitions, or resonant neutrino oscillations - but in this case 
resonances between coupled radio waves traveling in a tenuous medium. Interestingly, we find that sudden resonant 
transitions do nothing observable, while partially resonant transitions and the adiabatic limit produce very interesting 
results indeed. 

Application of electrodynamics to a magnetized neutral plasma (with an external magnetic field B = BB) gives 
£ ij = £ p$ij ~ iuB£ijkBk{d£p/ doj), with e p — 1 — uj 2 /oj 2 , where oj p is the plasma frequency, given approximately 
by up 1 = 4ire 2 n e /m e , and u>b is the cyclotron frequency uib — eB/m e . Small corrections of order kT/m e for a non- 
relativistic hot magnetized plasma || have no effect on the discussion. A common method of solution assumes that 
the magnetic field and the wave vector are parallel. One readily diagonalizes the dielectric tensor for this case, which 
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has transverse eigenstates of circular polarization E± and wavenumbers k + and fc_ for the same ui. It is tempting to 
focus on the part of B parallel to k, while treating other components (B±) as a perturbation. Then, goes the usual (but 
faulty) argument, the total rotation of polarization when the parameters vary along the line of propagation would be 
given by integrating (l/2)(fc+ — fc_) along the path, resulting in the astronomers' rule AO = RMA 2 for u> 3> u> p , [0- 
This seems rigorous enough in astrophysical circumstances, where us ~ 10MHz — 1GHz, n e ~ 10~ 6 cm~ 3 — 10 _3 cm -3 , 
and B w 10~ 9 G — 10 _5 G, and the approximation us 3> ujb,oj p is satisfied by several orders of magnitude. 

Our point of departure comes from noting that the propagation problem for B = is utterly degenerate, creating 
conditions for a corkscrew effect: the direction of B lifts the degeneracy and selects unique eigenstates. Treating the 
transverse parts of B as a perturbation will not do, because the effects of the perturbation are the same order as 
the splitting of eigenvalues. We solved the problem exactly for the two transverse eigenstates of Di = SijEj used 
in the adiabatic limit. One needs the eigenvalues of the inverse dielectric tensor in the 2x2 subspace orthogonal 

to k, which are given by e^ 1 = p^p — |g ""2 zf^ =t b ^ cos 2 ^ fc ( ? 2^2" The magnetic field direction has been 

represented by polar coordinates {dkB, 4>kB) in the coordinate system where k is along the z axis; b — lo^lob/^ <C 1 
is a dimensionless small parameter measuring the strength of the B field. For finite b, the system is non-degenerate 
for any direction of B. As cos(i9fcs) — > 0, the splitting of the eigenvalues reaches a non-zero minimum, exhibiting an 
avoided level crossing. 

The exact solution also reveals that the polarization eigenvectors deviate strongly from circular polarization when 
k and B are nearly perpendicular, or, more precisely, when cosies is of order b. A standard way to express the 
polarization of a radiation eigenvector uses the Poincare sphere jjj. In the hclicity basis, any normalized complex 
two-vector \s) can be associated with a 2 x 2 tracclcss, determinant —1, Hcrmitian matrix P, which satisfies P\s) = \s). 
If we write 

p_f cos?9p sin^pe - ^" 
I sm-d p e lct>p — cos-dp 

then the angles ($ p , (f> p ) define the polarization state of the vector. The positive and negative circular polarizations are 
represented by the top and bottom poles of the sphere; linear polarizations are represented by points on the equator. 
The solution for the eigenstate of definite frequency d\ is found to be 

tan(tfp) = 6sin 2 (?9 fe£ >)/cos($ fc£ ,); p = -2<f) kB , (3) 

with di following from orthogonality. To understand the physical effects, from (||), we see that the eigenvectors "follow" 
the magnetic field, but in a highly non-linear fashion. For b <C 1, and | cos(i9fc_B)| not too small, the eigenvectors are 
nearly circularly polarized, living on tiny circles with <j) p = —2<fikB near the top or bottom of the Poincare sphere. 
(The double-valuedness of this map is a reflection of the projective character of polarization). 

Something very dramatic happens when B rotates transverse to k. The dominance of circular polarization is upset, 
even for non-zero b arbitrarily small, as the components of B transverse to k receive all the responsibility for lifting 
the degeneracy. The position on the Poincare sphere jumps from its azimuthal position near the pole to the equator 
just as B z crosses zero. This will occur in any circumstance (see Fig. 1) when the sign of the component of the 
magnetic field in the z-direction changes (a " B z reversal"). 

It is convenient to consider the case where the magnetic field direction traverses a closed orbit along the path. 
We present the cases shown in Fig. 1(a), where the B rotates in a spiral about k without reversing, compared to 
another case [Fig. 1(b)] where B z reverses while undergoing an azimuthal twist. The corresponding trajectories 
of the polarization on the Poincare sphere are also shown. Although the system parameters return to their initial 
configuration after the cycle, the eigenstates pick up a geometrical phase, with magnitude on the positive (-E+) 
eigenvector given by the solid angle Q traversed on the Poincare sphere. The phase on the negative (E-) eigenvector 
is easily shown to be the same with opposite sign, so that the phase difference of the two circular polarizations is 
twice the solid angle swept out. This phase difference, which is just the angle by which the plane of polarization is 
rotated, is entirely of geometrical origin, and is in addition to the eikonal phase exp[i J kj(z) dz\. 

The geometrical character of the process can be intuitively visualized by considering the effects on a two-dimensional 
Cartesian coordinate frame (representing the set of polarizations) erected tangent to the sphere. As the coordinate 
frame is dragged along parallel to itself, its orientation accumulates information on the curvature of the state space. 
For the no B z reversal case [the spiral field of Fig. 1(a)], the solid angle fl is of order b 2 . A phase shift of order 
b 2 would be quite small, and give a rotation of polarization going like A 6 . On the other hand, suppose the field B z 
reverses, dragging the tangent frame from top to bottom of the sphere. At the bottom, the frame (oriented at the 
angle 4> p ) undergoes a twist by the azimuthal angle —2A4>kB, then B z reverses again, followed by +2A(f)kB- The net 
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phase shift is f2 = AAtfikB- The 4>b dependence of <f> p is crucial. Of course it is not necessary for b to be constant 
with z. Nor is a closed circuit necessary. When the slow parameters p k vary arbitarily, the corkscrew effect gives a 
rotation of polarization by an angle = 2 J p ^ dp k ■ A k (z)], which is of order unity and which is determined solely 
by changes of directions of B. Here A k = (1 — cos , dp)d k 4>p, where d k is the derivative with respect to the parameters 
[b(z), cosi?fcs(z), <fikB(z)], or, more directly, {dp,4> p ). (In Berry's convention [Q, the phases are already attached to 
eigenstates.) If there is more than one B z reversal along the trajectory, the various contributions track the history of 
the polarization states along the trajectory. 

Borrowing quantum terminology, we have just examined the adiabatic approach to resonance where a phase shift 
of geometrical origin occurs. Under adiabatic conditions, the "probability" (amplitude squared) of each eigenstate is 
conserved, while the phase-physics is non-trivial. The adiabatic approximation applies when pT 1 (z)dp(z) / 'dz <C kb. 
This means that parameters vary over a length L such that w^wbL/w 2 > 1, that is, that the wavelength is not too short. 
Comparing the scales, this is also the circumstances where Faraday rotation is non-negligible. Numerous interesting 
phenomena also occur outside of the adiabatic limit. In non-adiabatic resonance transitions the probability can hop 
from one state to another, again generating non- trivial dynamical phases. In the extreme case of the "sudden" 
approximation the parameters vary so rapidly (on the distance scale of the wavelength splitting) that the system 
jumps abruptly from old eigenstates to new eigenstates based only on the overlap (dot-product) of the states. For 
the magnetized plasma, this is the limit of co — > oo. With b (the dimensionless level splitting parameter) going like 
1/lu 3 , the splitting becomes so small that any physical variation of parameters becomes "sudden." Transitions occur 
on the basis of maximum overlap, and when the eigenstates switch, the probability switches correspondingly. This 
answers the paradox that the i?-field actually does lift the degeneracy of a gamma ray, but with negligible results: the 
corkscrew effect, to be sizable, occurs when the ordinary Faraday rotation is sizable. One cannot avoid the conclusion 
that under a wide variety of circumstances, a S z -reversal will produce non-trivial phase shifts between the helicity 
states and a substantial frequency-independent rotation of polarization. 

Let us now turn to astrophysical implications of these results. Current knowledge of intergalactic magnetic fields is 
scanty. Even the basic origin of established galactic fields is unknown. There are two major theories invoking internal 
dynamos or pre-existing primordial magnetic fields |J. Indications of a coherent field coming from excess RM's of 
extragalactic sources has been claimed, but the interpretation of the RM data is ambiguous ||. The existence of 
magnetic fields in galaxies at large redshifts has been argued to be a challenge to galaxy dynamo theory [Q. 

Recently, an indication of anisotropy in the propagation of radio waves on cosmological distances was reported 
||. After the Faraday effect was subtracted, a signed difference of angles called (3 was observed to behave as (3 = 
r/2A s cos(7) [In the analysis of ||, care was used to employ a variable (/?) which properly reflects the projective 
nature of polarization observables. Combinations such as x - ^ \x ~ ^1; an d the acute angle between \ an d ip cannot 
serve to detect cosmological birefringence.] Surprisingly, the direction of the anisotropy axis extracted in [^) coincides 
within errors with the direction of the dipole component of the cosmic microwave background Q , perhaps indicating 
that some type of intergalactic or cosmological medium may be at work. The possibility of systematic bias in the 
data was discussed, with emphasis on corrections made by the observers for our own galaxy. Studying the internal 
structure of the data, it was noticed Q that strong correlations were observed for about half the data for galaxies 
with z > 0.3, while no significant correlation was observed for the half of the data with z < 0.2 — 0.3. If something 
not taken into account were due to our own galaxy, it should affect all data for the extragalactic sources in the same 
way. As reported in Q , unexpected physics including the corkscrew effect that might occur in our local neighborhood 
cannot account for the observations. 

Whether a magnetic field of cosmological extent could explain the same correlation is a very interesting question. 
The turning-on of the correlation for z > 0.3 may be an indication that a system enormously larger than our local 
group trails a magnetic field or some other condensate of huge size. Galaxy density distribution studies tend to support 
the presence of coherent structure over vast scales, although this is controversial j[o|| . Complementary information 
might come from measurements of the cosmic microwave background polarization affected by primordial magnetic 
fields existing long before galaxy formation jy]]. Moreover, there could be a wide range of interpretations for the 
results of including anisotropic cosmologies of general relativity, an axion-like condensate domain wall, universal 
rotation, spontaneous rotational symmetry breaking, or antisymmetric degrees of freedom in gravitational theory [fL2|| . 
We cannot begin to answer all the questions here. A systematic approach using all the tools available, and orders of 
magnitude more data, would be able to make a significant scientific contribution. 

While focusing here on the parameters affecting linear polarization, there is clearly a scientific gold mine in the 
accumulation of all radio polarization observables which have never been interpreted taking into account the effects we 
discuss. It is very important to have data at different frequencies, in order to exploit the full potential of observations. 
The corkscrew effect might serve as a valuable discovery tool when larger data sets become available, with many 
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sources examined over a range of radio frequencies, supplemented by optical and ultraviolet polarization data. 
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FIG. 1. Cartoons illustrating a phase change from slow variations in the direction of B along the path of the wave, (a) B 
varies in a spiral, producing a tiny circular orbit on the Poincare sphere; (b) B z reverses while B$ twists, producing a flip and 
twist on the sphere. Corresponding points (1-5), with the configurations at 5 equal to 1, are shown. 
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